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ANNALS OF MATHEMATICS. 

Vol. IV. August, 1888. No. 4. 

ON THE DIFFUSION OF HEAT IN A HOMOGENEOUS RECTAN- 
GULAR MASS, WITH SPECIAL REFERENCE TO BARS 
USED AS STANDARDS OF LENGTH. 

By Mr. R. S. Woodward, Washington, D. C. 

I. Considering the vast amount of time and labor which have been spent in 
developing the science and art of comparing standards of length, and the won- 
derful precision at present attainable in such work, it would seem somewhat 
surprising that more progress has not been made toward a complete theory of 
the behavior of metallic bars under varying temperature conditions, [f we look 
over the scientific literature of the last half century, we shall find elaborate 
reports and memoirs on the appliances and methods adopted in determining 
coefficients of expansion and relative lengths of standards. But the chief object 
of all these appliances and methods has been to secure particular conditions, 
such as a steady low or high temperature of the bars under comparison; and the 
precise results for expansions and differences of length come out without refer- 
ence to the intermediate states and effects arising from the slow diffusion of heat 
through the bars while changing from one steady temperature to another. 
Indeed, a dexterous avoidance of all the more complex phenomena presented 
by a cooling or heating bar has been the characteristic feature of experimental 
success. Our authorities on standards of length give us multitudes of facts and 
figures concerning the history of a bar during particular periods, but the history 
begins and ends abruptly and we are left to conjecture concerning the facts and 
figures of the more interesting periods. 

The reason for this apparent anomaly, this paucity of information relative 
to the more recondite thermal properties of a bar, is obvious enough when we 
reflect that the present high grade of efficiency in the management of standards 
has only been attained after prolonged and diligent research. Obstacles quite 
independent of the standards themselves, of the most perplexing and baffling 
character, have had to be overcome. It is but natural and logical, therefore, that 
experimenters should have sought to avoid complexities as far as possible and 
direct their efforts so as to secure, in practical work particularly, the requisite 
precision at the minimum cost. 
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But precision is merely relative. What satisfies us to-day is almost sure to 
be unsatisfactory to-morrow. The perfection of methods and appliances which 
enabled us to overcome old difficulties has also revealed new ones, difficulties 
intimately related to, if not directly dependent on, the very states and conditions 
we have hitherto found it convenient to avoid or ignore. At present, precision 
of comparison appears to have reached a superior limit along the old lines of 
research ; and future progress seems possible only by an advance in our knowl- 
edge of the more intricate thermal properties of metals. Apparently we must 
recur to the extraordinary studies of Fourier and Poisson ; we must extend their 
work and adapt it to the practical needs of the computer ; and we must learn 
how to determine readily for any standard the constants which define its intrinsic 
thermal relations and the relations it sustains to media in which it may be placed. 

2. It is proposed to consider in the following pages the problem of a rectan- 
gular bar cooling or heating in a medium whose temperature remains sensibly 
constant. This is probably the easiest of the problems presented in the work of 
comparing standards of length, but it appears to be also the one whose solution 
will be most useful in a practical way. 

Much, if not most, of the preliminary work has been done by Fourier * and 
Poisson f, and the reader is referred to their treatises for all elementary details. 
In his Chapter VII Fourier has considered the case of a rectangular prism of 
infinite length, subjected to a constant temperature at one end, and exposed to a 
uniform current of air at zero temperature; and in his Chapter VIII he has given 
a pretty full treatment of the case of a cube cooling in air or any other medium 
of sensibly constant temperature. Poisson in his Chapter IX has discussed at 
great length the case of the distribution of heat in a bar whose transverse dimen- 
sions are very small, and in his Chapter XI he has briefly considered the case 
of a rectangular mass of unlimited dimensions. Neither of these writers, how- 
ever, appears to have had any practical applications in view, and elaborate and 
elegant as their work is, there remains much to be done to render it available for 
numerical applications. 

We shall assume that initially the bar has a uniform excess in temperature 
over that of the surrounding medium ; that the conductivity and thermal capacity 
of the bar and the emissivity are constant within the range of temperature con- 
sidered. 

Let the length, breadth, and thickness of the bar be denoted by 2a, 2b, and 
2c, respectively; and let any point within or on its bounding surface be defined 

*Thiorie Analytique de la Chaleur. Paris, 1 822. Analytical Theory of Heat, by Joseph Fourier, 
translated, with notes, by Alexander Freeman, Cambridge, 1878. Our references in the text are to Free- 
man's translation. 

\Thiorie Malhimatique de la Chaleur. Paris, 1835. 



IN A HOMOGENEOUS RECTANGULAR MASS. 103 

in position by a system of rectangular co-ordinates, whose origin is at the center 
of the bar and whose axes of x, y, z are parallel respectively to a, b, c. Call u the 
excess in temperature over that of the surrounding medium of any point x, y, z 
at any time t after the initial epoch. Then the partial differential equation, con- 
necting /, u, x,y, z, which obtains in all cases of the diffusion of heat, is* 

^ == *[a? + .^ + ^?J- (I) 

In this, k is the conductivity of the bar divided by its thermal capacity per unit 
of volume. 

In addition to the relation (i) pertaining to any element of the bar at any 
time (, there are three other equations which express the equalities between the 
heat conducted across any element of the faces of the bar and that carried away 
by the surrounding medium. If we call the conductivity K and the emissivity 
H, these relations are 

=P K-z- = Hu for x = ± a, (2) 

5u 
^ Kty = Hulor y=±b, (3) 

9u 
-t- K-7T- — Hu for z = ± c. (4) 

Finally, there is the initial relation expressing the fact that at the beginning 
of cooling all points of the bar have the same temperature. If u be the initial 
uniform excess in temperature of the bar above that of the surrounding medium, 
we have 

u = u for t=o. (5) 

The analytical solution of the problem consists in finding a function 

u = F(u ,x,y,z, t), 

which will satisfy all of the relations (1) . . . (5). 

3. The complete integral of equation (1 ), for the present case, as shown by 
Fourier in his Chapter VIII, consists of the product of three infinite series and 
the initial excess of temperature «„. We may express this integral thus: 





u=u V'V"V". 


»=00 

« = o 


»=00 

V" = I v n ", 
» = o 



V = I v n ', V" = 1' v n ", V" = 2' v,!"; (6) 

«=o »=o »=0 

*See Fourier's Chapter II, Theorem IV. 
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wherein 

vj = A H e - *'V cos /„*, 

v n " =B n e-^n'cos n„y, 
v n '"=C,e'- ktv n' cos v n z. 

In these expressions A,„ B n , C„, X n , //„, v„ are constants varying from term to 
term with the index n, which is zero or any positive integer. These constants 
must be such that (6) will satisfy (2), ... (5) as well as (1). 

For brevity, put k=H/K. Then, if we substitute u and its derivatives from 
(6) in (2), (3), and (4), there result 

aX a tan a/,„ = ah = ?, say, 

b/i n tan b/i n = bh = yj, say, (7) 

cv n tan cv n =ck = C, say. 

Since the values of c, y, and £ are known, the equations just written will deter- 
mine the constants X n , //„, v„. Each of these equations, it will be observed, has 
an infinite number of real roots, the successive values of which are to be used in 
the several terms of v n ', v n ", v n '" in the order defined by the index n. 

The constants A n , B m C n are determined in the following manner: Equation 
(5) requires that when / = o, 

n =00 « =00 « =00 

• I v,!= 1 V n "= I Vj"=l, 

tt = o n =0 « = o 

that is, considering v n ' alone, 

1 = A 9 cos l^x -\- A x cos XjX + A 2 cos X^x + • . . . 

If we multiply both sides of this equation by cosX n xdx and integrate 
between the limits o and a, it will be found that* 

A 4 sin aX n 

" 2aX n + sin 2aX n y -°> 

Corresponding expressions for B n and C n may be written down by an obvious 
interchange of symbols. 

4. It is seen, therefore, that the use of equations (6) as they stand, requires 
the roots (or a certain number of them) of the transcendental equations (7), these 
roots entering respectively into each of the three factors of the corresponding 

*See Fourier, Chapter VII, p. 313. 
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term of v „', v n ", or v n '". Herein lies the chief obstacle to the application of (6), 
and the principal feature of our work consists in showing how this obstacle may 
be overcome. 

We shall show first how the roots of (7) may be found, though our solution 
of the problem requires in any case only the smallest of these roots. 

For brevity, let us write 

<f,,=fyn, (9) 

Then (7) become 

0„ tan n = c, 

(.0) 

It will suffice to consider the first of these equations, since they are all of the 
same type. 

£z=ah is of zero dimension as respects length*, and it may, in general, have 
any value between o and 00 . 

In the case of an iron bar cooling in air, h may be about ^fft w 'th the cen- 
timeter as unit of length; so that for a bar of iron one meter long, we should 
have approximately (roughly, perhaps) 5 = ^. But for other media, h may be 
so large as to make ? exceed unity when a is no greater than a half meter. In 
any event, $ will exceed unity if a is sufficiently great. It becomes essential, 
therefore, to distinguish two cases; namely, that in which ?<i and that in 
which f >i. 

First, consider the case wherein £<i,with a view to expressing d n in a 
series of ascending powers off. When £ = oJ, we must have Q n = mz; i.e. n is 
either zero or a multiple of n. Differentiating the first of (10), there results 

(tan0 n + 0, lS ec 2 n )^|=i. (11) 

This gives, for £ = 0, 

-ng-= 00 when n = o, 
= (»^) _1 when n~y.i. 

* See Fourier, Chapter II, p. 1 28. 

fSir Win. Thomson gives (Encyclopedia Britannica, 9th Edition, Article Heat) H= 0.0002 for a 
polished copper surface, and A"=o.i6 for iron. Assuming the above value of #to apply for iron as well 
as copper we have k = jj ff . 

\ The physical meaning of this condition is that there can be no escape of heat from the bar. 
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It appears, therefore, that we cannot apply Maclaurin's series to the development 
of the first value of 9 n . But if we substitute for tan 9 n its value in series, we have 

f = 0. , + R 4 + -iW+ • • •; 

whence, by reversion, we get for the first roots of (10) 

9* = f - ic 2 + &? - &? + r 6) * 

<fo 2 = V - W + fa 3 -j& ft + T„ ( 1 2) 

h^t-tf + te'-tfeC+Ty 

These series will suffice for the calculation of the first values of 9 n , <p n , <p n . The 
convergence is rapid for small values of?, fj, Q, and quite adequate even when 
they are as great as unity. 

For the larger values of 9 n we may apply Maclaurin's series. Thus, by 
eliminating the trigonometric forms from (u), we get 

whence, by an application of the theorem of Leibnitz, we readily find that, 
for £ = 0, 

^ = +(«r) \ 

|!=-2M" J , (13) 

-jp = + 1 2(«s-)- 5 — 2(«n-)- 3 , 



Hence, we write 

ff. = * ff {i +.fH- J -c»- 4 +?[2(ro)- s -i«] + T t }. (14) 

This series converges rapidly, especially for large values of n; it shows, in fact, 
that the larger the root the more nearly is it equal to nn. 

Second, when ?>i, we may write the first of (10) in the form 

cot 0» = 9 n £ 1 , 

wherein c 1 = ? _l . (15) 

Now the first of these equations gives for £ t = o,t 

* We shall use throughout this paper T u to indicate that the next term of the series beyond those 
given is a term of the «th order. Thus the next term in the first of (12) is of the 5th order with respect 
to the variable f. 

f The physical meaning of this is that the heat is instantly dissipated on reaching the surface of the 
bar, or escapes "freely." 
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d n = l(2n+ i)n = &, say, 
dd n _ „ 
c£- ~*> 

-gi= + 29, (.6) 

Therefore, in accordance with Maclaurin's series, 

fl.=*[i-e- 1 +^-(i-** , )f- , + 7;]- (*7) 

Except for large values off, this series converges too slowly to be of much 
practical value. Our process, however, does not require its use. We give it only 
to facilitate a numerical verification of the equivalence of the original and trans- 
formed expressions for V, V", V '. See §7. 

5. It will appear in the sequel that, for small values off, rj, and £, the series 
of (6) converge rapidly toward their first terms with the lapse of time from the 
initial epoch. On this account, it is desirable to have a ready method of com- 
puting v„', v ", and v '". The factors in these terms are all in convenient form 
except A , B , and C , which are implicit functions of f, rj, and Q, respectively, as 
defined by (8), (9), and (10). They may be expressed directly in terms off, rj, Q 
in the following manner. Taking A„ as the type, (8) and (9) give 

4 sin g 



A 



2# +sin2# 



Substituting for sin 9 and sin 20 o their values in series, we find 

Replacing d 2 in this by its value in (12), there result 

B»=i+h-Th-?+T z , (18) 

c =i+K-Th: 2 +r* 

From (6) and (9), we have 

V =Age' kt (^^ cos O ~, 

V(j " = B e - *'<*•/*>' cos <p - d > ( 1 9) 

v '"= Cqe -*«*-")' cos ^ -• 



108 WOODWARD. ON THE DIFFUSION OF HEAT 

Computation of numerical values from these expressions is now easy, since 

<V ?o. <Po are s iven b y ( I2 ) and 4» B o> Q b y ( l8 )- 

In addition to v ', v ", and f '"-we need 

«=oo tz=oo n = oo 

I v,,', I v„", and I v n '". 
«=i «= i »= i 

These we proceed to expand as functions of £, rj , and £, respectively, by Mac- 
laurin's series. 

By reference to (8), (9), and (14) it is seen that when »>i and ? = 0, 

e — *«„* — ^ — kt (""/")* 

x 

cos A„x = cos tin- ; 

a 
d - d £?= + i2{-iy+\nn)-\ 



and -j£ = + (amt)-\ 



ds 



Hence, when £ — O, 



^ = — 2 (—'P 1 , -«(-/»)• cos«^f, 
a? jt w' a 

d*vj I2(— i) n+1 *««„/-« x 

+ 8 M- I )" +1 g -^(^/«)' cos «;r£ 
+ 4*(— ^0^ -m^^)' sin wr-, 



If now we write for brevity, 



C' = 4 2' ^ — H — ,-*«.«/«)• cos W7r f, 
<r' = 4 I' ^— ^ — *-*«»»/»>* cos «;r* 



(20) 
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we have for $ = o, 






9 


y vj 


n 


= I 




5c 


n 


= 00 


s 2 




n 


= 1 



4- ■> 
5P ^ 



^•'-'§'-f} 



Finally, if we denote by Si", S 2 ", etc., and by 5,'", S 2 "', etc., the series 
corresponding to (20) when a and x are replaced by b and _y, and c and £, 
respectively, there result 

V> = v > -S/f + ( 3 S/ -t^§- -^] f2 + r a. 
P' == V - S/'rj + [ 3S 2 " - I ^ ~*^f ] f + r s , (21) 

f'" = v>» - src + [ 3V - * ^f-' - * s -^p\ ? + r 3- 

Equations (12), and (18) to (21) afford a complete solution of the problem for 
the case wherein £, '<?, Q are each less than unity. The series in (20) converge 
with such rapidity that for any but the earliest stages of cooling their first terms 
will suffice. For practical applications to a standard of length of ordinary 
dimensions, cooling or heating in air, we may neglect all terms after the first 
in (21). 

The values of 5/, etc., in (20) might be expressed by a series of definite 
integrals similar to that derived in §8, but it does not appear that any practical 
advantage would be gained thereby. 

6. Considering still the case in which £, rj, and £ are each less than unity, it 
will be of interest to have expressions for certain average temperatures, namely: 
(a) the average temperature of the whole mass of the bar, (/3) the average tem- 
perature of any face of the bar, and (j) the average temperature along its longi- 
tudinal axis. 

(a). The average temperature of the whole mass is defined by the definite 
integral 

+a +i +c 

(8a6c) ~ l J udx \ dy\ dz. 
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A glance at (20) and (21) shows that this is equal to 

+a +i +c 

u {%abc) ~ ' J vjdx I v "dj J v '"dz + T 2 ; 

— a — b — c 

whence by reference to (19) we get the following result:* 

(») = ^, 5in v:y in ^ -^ + -• ^ + t. 

"o YD YO 

which, by reason of (12) and (18), is easily reduced to 

(a) = u^e ~ *«*»' ; - "»' + "- 2 > + T 2 . (22) 

(/9). The average temperature of any face of the bar, for example the face 
2a X 26, is expressed thus : 

- a ^b 

oJ"{4ab)- x judxCdy, 

— a -b 

wherein a'" stands for the second member of the last of equations (21) when 
s = c. Integrating, there results 

(p) = u oJ"{e ~ kt ^-' + *•*> + T 2 ). (23) 

It is easily seen that this differs from (22) by a term of the first order in £. 

(y). In a similar manner we find for the average temperature along the lon- 
gitudinal axis of the bar, 

(r)=«b *."*."'(<-**•" + 73, (24) 

in which <r " and <r '" stand respectively for the second and third of (21) when 
y = o and z = o. Obviously, (24) differs from (22) and (23) by terms of the first 
order in i] and £. 

7. Returning now to the case in which £, -q, and £ are each greater than 
unity, we proceed to show how the sums in (6) may be expanded in series of 
ascending powers of £ _1 , ij~ l , and £~\ Referring to (15) and (16) and using 
for brevity a, /?, y as defined below, we get for f _1 = f 1 =o, 

a = A„=+2(— i)"#-\ 

wr "■ 
£—<->*. 

» 

* Compare Fourier, Chapter VIII, p. 328, for the case in which a = b = c. 
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Ill 



„-*tt » = 



_f_^-*<<*/«)* ( 






= + 2kt{»lafe- u W, 



§* = — f6*/(#/«)»— (2^) 2 (^/«) 4 l^-*'(*/">', 



r = cos l,* = + cos #-> 






■h = -p 



(9 



cos «?- — 2«?- sin i?- > 
a « a 



If now we put PJ,' for P when £- l = f 1 = o, the preceding relations enable 
us to derive readily the following : 

K' = =2" i — r -^-*« 2 » + ">*<"-"«>' cos 4(2« + i>r-, 

JT 2»+I 2V ' ' « 

» = o 






$v 



&? 



W 



. 9VJ.9VJ ?V ' 



(25) 



Hence by Maclaurin's series 



SV 



V = V> +^f-'+ .... 



(26) 



in which V ", V '" and the derivatives of V", V" with respect to y~ l and £ _1 
are obtained from (25) by replacing a, x by b, y and <r, ^ respectively. 

The applicability of (26) for numerical computations will depend chiefly on 
the rapidity of convergence of the series expressing V ', V ", and V '". These 
series (see the first of (25) ) converge rapidly when the time is sufficiently great, 
but rather slowly during the earlier stages of cooling. They may be expressed, 
however, by a series of definite integrals which converge with extreme sud- 
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denness for all but large values of the time; and as the process of transformation 
from the former to the latter series is of great utility in other problems as well as 
in this, we shall give the steps of its present application rather fully in the next 
section. 

8. If in the equation * 

g-v/a.r^^L. C e —-'i'd r cos ^ 

o 

we put 

2 I 

and j3 =( 2 «+i)^> 

the value of V ' in (25) becomes 

F '= — 4-j3 Ce-y'^ t dr n i M {~ 1 J > , cos i(2» + 1) it- cos i(2» + 1) n X - 



2 r> * — °° (— iV" 



COS i(2«+ O^ [S*J 
COS |(2«+ i)tt f ? -^J 



•(27) 



Now the series under the sign I in this equation is periodically equal to zero 
and \tz, as y varies over its range from o to 00 ; and hence V ' consists of a defi- 
nite integral whose successive parts are multiplied alternately by o and \n. 
These facts may be shown in the following manner : 

Take the difference between the well known equation 

1— 2g cos p+g* = * +•£" C0S P+i i cos 2 P +Z 3 cos 3/ + • • • 
and the equation which follows from the substitution of n +p for/. The result is 

i [r=^ ¥7T? ~ \ + 2 g cos 7+p] =* cos / +^ 3 cos 3/ + • • • • 

For / in this, substitute go Q — (G +/). We find 

,r *-** i~g* I _ 

T Li -2^- sin (C +/) +J- 2 1 + 2 g sin (<?+/) +^J - 

g sm {G + p)-g 3 sm ?,{G + p) + . . .. 
*See Vol. Ill, p. 78, of the Annals of Mathematics. 
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Multiply this last by dt> and integrate between the limits — /"and n -\- P. The 
result is 

1 c r i — s* i — &* 1 

Z J Li -2g sin (G+p) + g* ~ 1 + v sin (£+/) + g i J *= S ' Say " (28) 
— /> 

But 2 is evidently identical with the quantity under the sign 2' in (27) when 

g= 1. 
S = ~' (29) 

p=—. 

2a 

The element-functions of the integrals in the first member of (28) are obvi- 
ously zero for g = 1 except when G -f-/ is ^(42 '+ i)t in the first and ^(42 : -\- 3) 7T 
in the second, i being zero or any integer. In these exceptional cases the ele- 
ment-functions are infinite; and the values of the integrals depend on the limits 
— Pand P-\-it only so far as they determine within what range y can give to 
G -j-p the above named multiples of \it. Hence, in evaluating either of the 
integrals, as for example the first, for any case in which sin (G -\-p)= -\- 1, we 
may enlarge the limits to — 00 and + 00 without affecting the result. This 
premised, the required evaluation may be accomplished by a process frequently 
employed by Poisson in the treatment of similar integrals.* 

Substitute for g, 1 — w, where w is a small quantity approximating zero as 
g approaches unity. Also, replace G -\- p by ^(41 '■+ i)n + 1 in the first integral 
and by £(42 + 3)^ + q in the second, q being, like w, a small quantity. Then, 
neglecting terms of the third and higher orders in q, either of the integrals 
becomes 

+ t, 
! r wdq 

wherein q l and q 2 , as seen above, may have any values greater than an infini- 
tesimal. Hence the value of either integral is 

A (arc tan ±-2l_ arc tan J=1 ^ i 1 =in 
when w = o. 

*See, for instance, Chapters VII and VIII of the Thiorie Mathlmatiqut 4e la Chalevr. 
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We must now find the limits of G, or rather y as defined by the second of 
(29), within which the factor Q is \n. For the first integral the lower limits of 
y are given by the equation 

G = U4i+i)"-n ~P 
or y = (4/+ \)a — 2a — x; 

whence 

y = — ( a-\-x) for i=o, 

= + ($<* — *) for«'=i, 
= + (ja — x) for i = 2, 



and the upper limits are given by the equation 

whence 

?-= a-\-x forz'=o, 
=:c,a + x for z'= 1, 
= ga + x for 1 = 2, 



Similarly, the lower and upper limits for the second integral are, respectively, 

a — x and 5a + x for z'= o, 
5« — .r and ja-^x for z = 1 , 



Hence, bearing in mind that j- is positive only, we find for the transformed 
expression of (27) 

a-\-x 3a -f x $a-\-x 

o m, — x -\a — x 

(30) 

a — x 3a — x 5a — x 

o a+x y+x 

If now we write for brevity 

r„ = a/2\/(kt), 
r =x/2\/{kt), 

and replace y by 2yy{ki), there results 



(30 



v >' = y-lf e - y ' d r-f e - y ° d r +••■]■ (32) 
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The corresponding expressions for V " and V '" are obtained by merely inter- 
changing a, x with b, y and c, z. 

It will be observed, of course, that the series of definite integrals in (30) or 
(32) possesses all the properties and requirements of the less complex series in 
the first of (25). We shall not trace out these details, but the reader's attention 
may be drawn to the remarkable way in which (30) and (32) become zero for 
x= ±a and unity for t=o. 

In case £, rj, £ are so large that we may neglect all terms after the first 
in (26), 

u = u,V,'V a "V,'". 

This is always an approximate expression for the temperature when £ , rj, 
and £ are greater than unity ; and, other relations remaining constant, its degree 
of approximation increases as the dimensions a, b, c increase. For points near 
the surface of a very large mass, comparable in dimensions and other circum- 
stances with the earth, say, we may write 

(a-x)/*y\kt) (b-y)/iy/(kf) (c-z}/i^(kt) 

u = 8«„(ff)-» C e-y°d r f e-y'dy f e—*'dy. (33) 

000 

For points near the surface of a slab of infinite dimensions in the directions 
of b and c, and finite but great dimension in the direction a, the last equation 
gives 

(a-x)/^(kt) 



u=z^s. Ce-y'dr, 



a well known result. 

The derivatives of V a ' with respect to x and t which enter the last two of 
equations (25), are easily found from (32). For brevity, write 

/;, = 2tir + r — r, 
p 2 = 2nr + r + r, 

Hi = K 1 "(— i)» [>-*■ — e-p*°l 
n = o 

R 2 = n Z X (- i)« [p x e->s + P2 e-^l ( 34 ) 

« = o 

* 3 = *<M- i)"L>i , '-"'-ft , '-"1. 



n = o 

n = 00 



* 4 = 1 (-i)-[> 1 v-«'.- + ft ^-p.']. 

n = o 
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Then we have 

W _ R* _ 

dx 2 kty'iz 

9V '___ R 2 



(35) 



3t ty'lZ 

dxdt ' ^ 2tv\nkt) ^ ~~ 2K *>- 

9. We are prepared now to derive any species of average value of the tem- 
perature for the case under consideration. Of these, two may suffice, namely : 
(a) the average temperature of the whole mass and (/J) the average temperature 
wer any face of the mass as 2a X 2b. 

(a). For brevity put 



M' 



7 ^{2a)- l jVdx > 



+ 6 



Then 

If now we write 



M" =(2b)- l $V"dy, 
M"' = {2c)-^V'"dz. 



(a) = u a M'M"M'". (36) 

+ a 

N = (2a)- l fv 'dx, 

~"+« (37) 

— a 

it appears from (25) and (26) that 

M' = N+NJ~ l +T 2 . (38) 

Similar expressions obtain, of course, for M " and M'". 
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It remains to evaluate the integrals in (37). They are simple in form but 
present some of the difficulties peculiar to the class of functions to which they 
belong. The first of them is readily evaluated by direct operation on the first 
of (25), while the second comes by a more tedious process through the inter- 
vention of the expressions in (31), (34), and (35). The results are 



JV=8tt- 2 I 



e — it(?n + i,*ir*/ta* ) 



(2«+i)r 



„ = o (2» + O 2 

N 1 = Tt-* n ~s'°(— i)»rr - 1 (^-t» v o'_^-4(« + i)v o ») + 2 (2« + i)§e-?dp\. 
The steps in the derivation of N x are briefly these: 

+ a +a 

( 2a ) - 1 r£^Zj dx = -L- CCR, + —4r\ R&) dx 

— a — a 



(39) 



r. f. 



2(0 + i)r„ 



= r B - l 7r-i I [R3+ 2 I (- i)» fe-»'(2nr + r - ^ 

I -r. «— o J 



I 



72 = 00 

= r -V-» 2' (-i)« 
» = o 



2(» + i)r„ "\ 

[*-<^+2(2»-|-i)r {*-^ 



whence the expression in the text. 

The second members of (39) though complex in form are very easily evalu- 
ated, by reason of the rapid convergence of the series on which they depend. 
It maybe noted that for t=o, N=\ and iV 1 = o,and that for t=<x> , N== 
7V, = o, as required by the conditions of the problem. 

(/?). To get the average temperature of the face 2a X 2b, we have only to 
multiply u M'M" by the value of V" in (26) when z = c. Calling this value 
V e '" we have 

(fi) = u ly M'M"V c '". (40) 

10. It is believed that the solutions given by (6), (21), and (26) will answer 
most of the requirements in practical applications to standards of length of rect- 
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angular form. It is possible, however, that the developments with respect to 
$ , ij, and £ may prove to be in some cases insufficiently convergent. This can 
happen only when $, y, and £ are near unity, and the obvious modification requi- 
site for the treatment of such cases, is a development with respect to I — $, I — 7], 
and I — £. But in the absence of a definite knowledge of the circumstances 
giving rise to this phase of the problem, it does not appear useful to do more 
than suggest the modification. 

1 1. For the purpose of illustrating the application of the preceding theory, 
we shall work out in detail some numerical examples. 

First, take the case of an iron bar, \ m — ioo"" long and 4""' X 4"" in cross 
section; cooling in air. We have for iron, according to Sir W. Thomson,* in 
C. G. S. units, £ = 0.185, and we may, as in §4, assume k = -g^. Then, since 
a — so"™, b — 2™, and c = 2 cm , 

£ = tV» r i ~ C — T5T- 
With these data, equations (12) give 

log O 2 = 8.78689 — 10, 
log fo 2 - 7.39758— 10, 
log$V=-- 7-39758— 10. 



Also 



and 



From (18) we find 



d = 14° 10' 35", 
fo= 2 51 49, 
^0= 2 51 49; 

log k \ — \ = 4.65621 — 10, 
log k £b = 6.06269 — IO > 
log k \ — \ = 6.06269 — 10. 

log A = 0.00440, 
log i? = 0.00018, 
log C = 0.00018. 



Encyclopedia Britannica, 9th Ed., article Heat. 
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We are now prepared to compute the factors V, V", and V" of (21) for 
any values of t, x, y, z. For example, let us determine these factors and hence 
the temperatures for several points along the axis of the bar at the end of twenty 
minutes from the initial epoch. That is, make /= 1200"; x successively O, \a, 
\a, \a, and a; and y — z = o. From (19) and (20)- we get by means of the 
preceding data the following table of values for use in (21): 



X 

a 
O 


w 


s/ 


3^' 


as,' 

1 St 


3SJ 

9x 


1 .0047 


+ 0.0828 


+ 0.0255 


+ 0.0075 


+ 0.0000 


\ 


1.0028 


+ 0.0596 


+ 0.0181 


+ 0.0053 


+ 0.0170 


\ 


0.9970 


+ 0.0015 


+ 0.000 1 


+ 0.0000 


+ 0.0268 


3 

3" 


0.9875 


— 0.0596 


— 0.0181 


— 0.0053 


+ 0.0208 


I 


0.9741 


— 0.0859 


— 0.0257 


— 0.0075 


+ 0.0000 



In addition we have 



and 



v " --= v '" ------ 0.87054, 

S," = S,'" = 4- 0.0828, 



I 









+ 0.0330. 



Hence we find the following table of values of V, in which the terms in f and f 2 
are v given separately in order to show their magnitudes relatively to v ' : 



X 

a 


■-1 < 

<0 


c 3? 


** IP" 


V 


O 


+ I -0047 


— 0.0052 


+ 0.00013 


= 0.9996 


1 
4 


4- 1.0028 


— 0.0037 


+ 0.00016 


= 0-9993 


i 


+ 0.9970 


— O.OOOI 


+ O.OOOIO 


= 0.9970 


I 


+ 0.9875 


+ 0.0037 


— 0.0000 1 


= 0.9912 


1 


+ 0.9741 


+ 0.0054 


— 0.00013 


= 0-9794 
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We find also that 

V" = V" = 0.87054 — 0.000207 + 0.0000002 
= 0.87033. 

It is now only necessary to form the products V V" V" in order to get the 
temperature u in terms of u . The results are 



t = 


2o m , y = z = 0. 


X 

a 




u 

»0 


0.7571 


1 


0.7570 


1 

2 


0-7553 


1 


0.7509 


I 


0.7420 



It appears from the values in the second paragraph above that the terms in 
f 2 , rf, and £ 2 are, for the assumed values of those ratios, already practically insig- 
nificant at the end of twenty minutes from the initial epoch. Hence in computing 
temperatures for greater values of the time we may properly neglect all terms of 
the second order. For the purpose of showing the march of the temperature 
with the time, we have computed the values which u has at the above named 
points along the axis of the bar for intervals of twenty minutes up to t= 5* 40"*. 
Since minute precision is not needed in this work, we have used V = v ' — £5/ 
and V" = V" =v ". The resulting values are shown in the table below. The 
average temperature of the whole mass of the bar obtained from (22) is also 
given in the last column of the table. This average is only slightly greater 
(tsW^ P art ) tnan tne average temperature of the face 2a X 2b* 

♦The question whether under similar circumstances there can be any detrimental internal strain due 
to inequalities in distribution of temperature in standards of length of ordinary forms and dimensions, 
appears to receive a decisively negative answer in the computed values of the table. Such standards are 
generally smaller in cross section than our assumed bar, and would, therefore, other things being equal, 
sooner attain a sensibly uniform temperature. Hence we infer that, so far as flexure due to temperature is 
concerned, it is immaterial where the graduations on a standard are placed — whether on its external surface 
or on its neutral plane. 
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Time 


Values of temperature u on axis 


of bar where 


xja = 


Average 

for 
mass of 

bar. 





1 
f 


\ 


t 


1 


h m 














o 20 


O757*o 


0757*0 


0.755*0 


O75l*o 


O-742« 


0754*0 


40 


0-573 


0.572 


0.570 


O.565 


0.558 


O.568 


I CO 


433 


0.432 


0.430 


O.426 


0.420 


O.428 


I 20 


0.326 


0.326 


0.324 


O.32I 


0.317 


O.323 


I 40 


0.246 


0.246 


0.244 


O.242 


0.239 


0.243 


2 OO 


0.185 


0.185 


0.184 


O.182 


0.180 


O.183 


2 20 


0.140 


0.140 


0.139 


O.I37 


0.136 


O.I38 


2 40 


0.105 


0.105 


0.104 


O.IO3 


0.102 


O.IO4 


3 °° 


0.079 


0.079 


0.079 


O.O78 


0.077 


O.O79 


3 20 


0.060 


0.060 


0.059 


0.059 


0.058 


0.059 


3 40 


0.045 


0.045 


0.045 


O.O44 


0.044 


0.045 


4 00 


0.034 


0.034 


0.034 


O.O33 


0.033 


O.O34 


4 20 


0.026 


0.026 


0.025 


0.025 


0.025 


0.025 

J 


4 40 


0.019 


0.019 


0.019 


O.OI9 


0.019 


O.OI9 


5 00 


0.015 


0.014 


0.014 


O.OI4 


0.014 


O.OI4 


5 20 


O.OII 


O.OII 


O.OII 


O.OII 


O.OII 


O.OII 


5 40 


0.008 


0.008 


0.008 


0.008 


0.008 


O.OO8 



12. Turning now to the question of numerical applications for the case in 
which $, yj, and £ exceed unity, it will be advantageous to first set down some 
additional formulas for computing the differential coefficients which appear in 
the second members of (25) and (26). 

To avoid writing frequently recurring complex expressions, let us make the 
following abbreviations : 

n --■■■ 00 
X l = 1 (— i)"(2«+ i)*>-("/-4'-.> a (*« + «>* cos£(2«+ 1) na~ x x, 
» = o 
w = oo 

X 2 = 2' (— 1)" e -<.*/*'*>• (" + 1 )' sin \{in + 1) na~ l x y . 
» = o 



X 3 = E (— 1 ) n {in + 1 fe - <"/•"■•>' ( 2 » + i )'cos{(2k+i)to- 1 x, 
n = o 

X±= 1 {— 1 )" {in + 1 ) 2 e - <*/•"•<•>' < 2 » + -) 3 sin \{in + 1 ) za - ' x. 



(41) 



» = o 



r *=a'/4kt. 
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Then, the second and third of (25) become 

|^ = + ^o" 2 X x + 2a- l xX 2 , 

W _ ff(l + fl _ v + | ro - i)Xi + ^^ -^ 3 <4 ' j 

— 4« - "^ + i Af 2 « ~~ ^- 

The equations just given will be useful when r is small, since the functions 
in (41) will then converge very rapidly. For the converse case, we must make 
use of the differential coefficients given by (35). Introducing these in the second 
and third of (25), we get 

??= + £<*+*> 

= + ??2. [y-^o -»■)*+ ^-('0 + r)>_2^-( 3 r„ -r)'_ . . . ], (43) 

j/TT 

^ = - ^ [W + O*. + »*. + *J- 

The last expression assumes very easily managed forms for certain values of r 
used in the sequel. For the benefit of the reader who may desire to verify our 
work, these forms are given below : 

™= — 16 l<L-(e-r°' —i$ e -9rs _|_ 65* -w>* — ...)forr = o, 

= — 4 -St- (y-«r„)», + 3 ^-9(Jr<,)»_ 25 ^-»5(Jr„)' _ _ _ ,)f orr = l r (44) 
>- 3 

= + 64 — °- (*— » r <>* — 4e- l6r °* + q^-sS'-o 3 — . . .)forr = r . 

The series in (43) and (44) converge with extreme suddenness for large 
values of r , and they cease to be applicable only when r is very small. 

With preliminaries thus arranged, let us consider the behavior of a bar of 
the same dimensions and material as the bar of the previous example, cooling in 
a medium which produces an emissivity great enough to make ?, tj, and £ exceed 
unity. Without assigning to these symbols any particular values, it appears that 
the bar would very quickly assume, the temperature of the medium. For we 
have, as in § 11, 

a = 50°", b = c=2 om , 

and £ = 0.185; 
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and since for the transverse section of the bar 

... P i 
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4kt kt 



the quantities V ", V a '" and the derivatives of V" and V" relatively to )? _1 and 
£ _1 become practically nil soon after the initial epoch. Thus, from (25) and (42), 
we find the following table of values which show the rapid decrease of those 
quantities with the time 



/ 


V" 

* 


dV" 


cPV" 
9r)-* 


8 

30 
60 
90 


O.O4I 

O.OOI3 

O.OOOO4 


+ 0.284 
+ 0.018 

+ 0.0009 


+ 0.991 
+ 0.195 
+ 0.016 



Hence we see that, with such values of £ , rj, and £ as render equations (26) appli- 
cable, the bar in question would have sensibly the same temperature as the 
surrounding medium in less than two minutes after the initial epoch. 

Other conditions remaining constant, the rate of cooling of a mass dimin- 
ishes very rapidly with an increase of its dimensions. In order to strikingly 
illustrate this fact by means of the analysis, we may contrast the cooling of the 
bar just considered with the cooling of a cubic meter of the same material. The 
use of the formulas will be sufficiently indicated by computing the temperatures, 
corresponding to a few values of the time, at the five points designated by the 
following values of their co-ordinates: 



X 


y 


z 



±\a 

± a 













To compute V, V" , V" for such values of the time as we have assumed, 
one may use equations (25) and (42) or equations (32) and (43), the latter being 
best for the smaller and the former best for the greater times. We have used 
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both sets of formulas as a check on the computed results. The quantities needed 
in equations (26) are given in the following table. V B ', V ", and V '" and the 
corresponding derivatives of V, V", and V" are equal for x=y = z = o. 
Hence the values in the first section of the table are those required by the 
second and third of (26). 



t 


X = 


0, y -— 0, z 


= 0. 


x= ?z%a, y = o, 


^ = o. 


K' 


9V 




V 


dv 


3*V 


h m 
O 20 


0.965 


+ 0.227 


-2.55 


0.764 


+ 0.940 


— 5-33 


I OO 


0.659 


+ 0.852 


— 3-50 


0.467 


+ 0.874 


-1.85 


2 OO 


0.342 


+ 0.905 


— 413 


0.242 


-f 0.825 


— 1.1 1 


3 00 


0.177 


+ 0.699 


+ 1.26 


0.125 


4- 0.587 


+ 3.80 






/ 


X 


= ± a, y — 0, z = 0. 






K 


/ 


9V 


3"V 


h m 
O 20 







+ 1.893 


+ 0.00 






I OO 







4- 1.042 


+ 0.84 










2 OO 







+ 0.537 


+ S45 










3 00 







-f 0.278 


+ 9-37 







Let us now suppose that f — y = £ = 10; z. <?., «^ = 50A = io, or ^ = 
H/K=\. Then we have for the centre of the cube when /= 20 ro , 

P = F' = V" = 0.965 + 0.023 — 0.0 1 3 + T g 

= o.975 > 
and u = (o.975) 3 « = 0.927 u . 

For a point on either axis of the cube half way between the centre and a face, we 
have when t= i\ 

V = 0.467 + 0.087 — 0.009 + ^3- 
V" = V" = 0.659 + 0.085 — 0.018 + r 3 ; 



and hence 
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u = 0.545 X (p.726)X — 0.287 «„. 
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In a similar manner the temperature at the centre of a face of the cube is 
found. The following table gives all of the temperatures derivable from the data 
in the preceding table. The values of the co-ordinates may, of course, be inter- 
changed, so that the temperatures given in the last two columns are common 
respectively to six points. 



Time. 


Temperature. 


/ 


At centre of 
cube. 


At middle 

point of 
semi axis. 


At centre of 
face. 


h in. 

20 

1 00 

2 00 

3 00 


0.927*0 
0.384 
0.069 
0.016 


0-79°*o 
0.288 

0.054 

0.013 


o.i8o« 

O.O57 
O.OI4 
O.OO5 



As a final example in this section we may give the average values of the 
temperature of the whole mass and of any face of the cube just considered. 
These values are easily computed from equations (36) to (40). Using the same 
values of the time as those in the last table above, we arrive at the following 
results, exhibited also in tabular form. The quantities A^ and N x are computed 
from ($g),M' from (38), and VJ" comes from the data in the third section of the 
second table above. It is assumed that $ = rj = £ = 10. 













Average temperature 


t 


N 


N x 


M 


VJ" 


















Of whole mass. 


Of any face. 


h m 

20 


0.689 


+ i-34 


0.823 


0.189 


o.56« 


O.X3*o 


1 00 


0.422 


+ 1-53 


0-57S 


0.113 


0.19 


O.O4 


2 00 


o.33i 


+ i-37 


0.468 


0.108 


O.IO 


0.02 


3 00 


0.1 13 


+ 0.99 


0.212 


0.122 


0.0 1 


O.O I 
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13. It appears desirable to indicate in the last section of this paper in what 
ways the conditions of the problem considered and the analysis derived may 
be useful in determining the thermal constants of a bar; namely, the thermal 
capacity of C, the conductivity K, and the emissivity H. 

The principal conditions of the problem we have treated are first, that the 
bar has initially a uniform temperature, and second, that it cools or heats in a 
medium whose temperature remains sensibly constant. The facts most readily 
observed under these conditions are the initial excess in temperature of the bar 
above that of the medium and the variation of its temperature or some function 
thereof with the time. From a sufficient number of such facts it would appear 
that all of the constants might be determined. To fix the ideas, suppose the 
excess in temperature at the centre of the bar or the average excess of its entire 
mass, can be measured at any time. Call this excess u'. Then we have 

u'=f(C, H, K); 
and if we write 

C a + AC = C, 

H a + AH=H t 
K, + AK=K, 

in which C , // , and K a are approximate values of C, H, and K found in any 
manner, we shall have the following observation-equation : 

^AC+^AH+^AK +/(C , H , K ) - u> = An'. (45) 

From a series of such equations, differing from each other sufficiently, one 
could by the method of least squares derive precise corrections to the approxi- 
mate values C , U a , and K . The applicability of this process, however, is limited 
narrowly by the difficulty of securing a sufficient range in temperature in a given 
bar, or by the difficulty of securing a sufficient range of dimensions in different 
bars or masses of a given material. 

Although it does not appear practicable under the special conditions of the 
present problem to derive all the thermal constants of a bar, there exist between 
those constants under certain circumstances some relations which may be accu- 
rately determined. These we proceed to point out. 

Firstly, when the emissivity and the dimensions of the bar are small, we 
have from (9) and (22), the average temperature of the bar's mass 

u' = u e ' *«> k °° + *»' + »•■)+ T it 
whence k (V + ft* + v 2 ) log e = t~ x log -» 
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But 

and from (9) and (12) 

V + ft* + V = (*-"+ 6~ l + c->) % 4- T t . 
Hence, to terms of the first order in H we have 

This expression gives the emissivity // in terms of the thermal capacity C when 
the second member is known by observation. Under the assumed condition- 
that H is small, u' may be safely inferred from the length of the bar (supposing 
it to be of small cross section). In case there are several observed values' of u' 
and the corresponding time /, let W be the most probable value of the second 
member of (46). Then, since u' or log u' only can properly be considered sub- 
ject to error of observation, we have in the usual notation of least squares 



W-. 



[< log J] 



From this and (46) therefore, 



(47) 



^ CK >-.+^-,-. (48) 

Secondly, when the emissivity is large or when £, rj and £ are such, we 
have from (6) and (26), 

£ = V'Vi'V'" [1 + [a-Wr* |S + • • •) § + T,J (49) 

Now V ', V ", V '" and the derivatives of V , V", V" involve the coefficient of 
diffusion or C and AT only; and hence when these are known an observation 
of u and t will give //through the relation (49). On the other hand, iff, y, and 
£ are so large that the terms beyond the first in (49) may be neglected, an obser- 
vation of u and t will give the coefficient of diffusion k. The method of treating 
a series of observed values of u and t to determine H or k, is sufficiently indicated 
by the remarks pertaining to equations (45) and (47). 



